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Abstract 

This paper describes a determination of the Cabibbo-Kobayashi-Maskawa matrix element \Vcb\ 
using the decay f£. We perform a combined measurement of this quantity and of the 

form factors , -Ri(l), and -^2(1) which fully characterize this decay in the framework of heavy- 
quark effective theory, based on 140 fb^^ of Belle data collected near the T(4S') resonance. The 
results, based on about 69,000 reconstructed B^ D*-i+ue decays, are = 1.293 ±0.045 ±0.029, 
Ri{l) = 1.495 ±0.050 ±0.062, i?2(l) = 0.844 ± 0.034 ± 0.019 and J^(l)|Kfe| = 34.4 ± 0.2 ± 1.0. The 
B^ D*^i^U£ branching fraction is found to be (4.42 ± 0.03 ± 0.25)%. For all these numbers, the 
first error is the statistical and the second is the systematic uncertainty. All results are preliminary. 
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FIG. 1: Quark-level Feynman diagram for the decay — > D* I'^v^. 

I. INTRODUCTION 

The study of the decay D*^i'^i'£ is an important item on the B physics agenda 

for many reasons. First, the total rate is proportional to the magnitude of the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element Vcb [?] squared. We can thus determine this 
quantity from the measurement of this decay. Second, B^ — > D*~£~^i'£ is a major background 
for charmless semileptonic B decays or semileptonic B decays with large missing energy. 
A precise knowledge of the form factors of this decay will thus help reducing systematic 
uncertainties in these analyses. 

Our analysis technique follows closely previous studies of this decay using e+e~ 
T{AS) data 0,311, i.e., we reconstruct the D* meson and the charged lepton only, without 
making any requirement on the other B meson in the event. The main difference to ear- 
lier analyses [2I, S] is that we measure the CKM matrix element and all three HQET form 
factors of this decay simultaneously. Also, by using a novel reconstruction technique of the 
B meson 4-momentum, we achieve a better resolution in the kinematic variables describing 
the B^ — i> D*~i^i>e decay, which translates into an improved determination of the HQET 
form factors. 



II. THEORETICAL FRAMEWORK 



A. Kinematic variables 

The decay B^ D*~i^i'e ^ proceeds chiefly through the tree-level transition shown in 
Fig. [H Its kinematics can be fully characterized by four variables: 
The first one is w, defined by 

^._ Pb- Pd* + ml. - 

where rriB and rriD* are the masses of the B and the D* mesons (5.2794 GeV and 2.010 GeV, 
respectively ^]), pb and po* are their four-momenta, and = {pe + PuY ■ In the B rest 
frame, approximately equal to the T(4S') center-of-mass (cm.) frame, the expression for w 
reduces to the Lorentz boost 7£). = E£)*/m£)*. The ranges of w and are restricted by the 
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FIG. 2: Definition of the angles Oe, Oy and x for the decay 5° ^ D*-£+iyi, D*- D^tt-. 
kinematics of the decay, with = corresponding to 

w^ax = ~ 1.504 , 2 

and Wmin = 1 to 

?max = i^B - mo- f = 10.69 GcV^ . (3) 

The point w = 1 is also refered to as zero recoiL 

The remaining three variables are the angles shown in Fig. [2) 

• di, the angle between the direction of the lepton in the virtual W rest frame and the 
direction of the W in the B rest frame; 

• 6y, the angle between the direction of the D meson in the D* rest frame and the 
direction of the D* meson in the B rest frame; 

• X) the angle between the D* and W decay planes in the B rest frame. 



B. Four-dimensional decay distribution 

The Lorentz structure of the B^ D*^i^i'( 
of three helicity amplitudes H_, and Hq), 



decay amplitude can be expressed in terms 
which correspond to the three polarization 
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states of the D* , two transverse and one longitudinal. For low-mass leptons (electrons and 
muons), these amplitudes are expressed in terms of the three functions /lyij(w), and 

R2{w) a 

ms „ r, J / hA,iw)Hi{w) , (4) 



where 



Hiw) 



Ho 



2Vl -2wr + r^ 
VI - 2wr + r2 (^1 ± 



w 



\Ri{w] 



1 



1 — r 

(W- l)(l-i?2H) 

1 — r 



(5) 
(6) 



with R* = {2y/mBmD^) / {rriB + niD*) and r = mD*/mB- The functions Ri{w) and R2{w) 
are defined in terms of the axial and vector form factors as, 



V{w) 



R2(w) 



W + 1 

Ri(w) 2 



Ai{w) 



By convention, the function (w) is defined as 

1 2 



AH 



AH 



(7) 
(8) 

(9) 



+ 1 

For w ^ 1, the axial form factor Ai{w) dominates, and in the limit of infinite b- and c-quark 
masses, a single form factor describes the decay, the so-called Isgur-Wise function jl, |9|. 
The fully differential decay rate in terms of the three helicity amplitudes is [lo| 



dw;d(cos6'£)d(cos6'y)dx 8(47r) 



- 1(1 - 2wr + r^)Gl\V, 



cb\ 



(10) 



X {{1 - cosOef sin^ OvHliw) + {1 + cos Oef sin^ By Hl{w) 

+ 4 sin^ Be cos^ OyH^iw) - 2 sin^ 6^ sin^ By cos2xH+{w)H_{w) 
— 4sin6'£(l — cos6'£) sin6'v' cos6'v' cosxH+{w)Hq{w) 

+ 4sin6'£(l + cos6'^) sin6'y cos 6'y cos xH^{w)Ho{w)^ , 

with Gf = (1.16637 ± 0.00001) x 10"^ GeV~^ By integrating this decay rate over all but 
one of the four variables, w, cosO^, cos 6'y, or x, we obtain the four one-dimensional decay 
distributions from which we will extract the form factors. The differential decay rate as a 
function of w is 

dr 



dw 48tt^ 



F 

m 



where 



(w + 1 

1 



D*{^B - rriD*) g{w)T'^{w)\Vch\'^ , 

1 — 2wr + 



(11) 



(1 -r)2 



X 



l + Riiw) 



w + 1 



, w — 1 



-\ 2 
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and Q{w) is a known phase space factor, 



g{w) = Vw^ - l{w + l f 



1 + 4- 



w 



1 — 2wr + 



w + 1 1 



In the infinite quark-mass hmit, the heavy quark symmetry (HQS) predicts = 1. 

Corrections to this hmit have been calculated in lattice QCD. A calculation, performed 
in the quenched approximation, predicts (including a QED correction of 0.7%) ^(1) = 



+0.030 
-0.035 



This value is compatible with estimates based on non- lattice methods |l2| 



0.919 

A recent unquenched lattice result, J-'{1) = 0.930 ± 0.022, is still preliminary [13 



C. Form factor parameterization 

The heavy quark effective theory (HQET) allows to obtain a parameterization of these 
form-factors. Perfect heavy quark symmetry implies that Ri{w) = R2{w) = 1, i.e., the form 
factors A2 and V are identical for all values of w and differ from Ai only by a simple kinematic 
factor. Corrections to this approximation have been calculated in powers of AqcD/'^b and 
the strong coupling constant a^. Various parameterizations in powers of {w — 1) have been 
proposed. Among the different predictions relating the coefficients of the higher order terms 
to the linear term, we adopt the following expressions derived by Caprini, Lellouch and 
Neubert Q, 

hA, (w) = hA, (1) [1 - 8p^z + (53p^ - 15)^^ - (231p^ - 91)/] , (12) 
Ri{w) = i?i(l) - 0.12(w - 1) + 0.05(w - 1)^ , (13) 
R^^w) = i?2(l) + 0.11(w - 1) - 0.06(w - 1)^ , (14) 

where z = {^/w + 1 — v^)/ {s/w + 1 + \/2). The three parameters p^, i?i(l), and i?2(l), 
cannot be calculated; they must be extracted from data. 



III. EXPERIMENTAL PROCEDURE 
A. Data sample and event selection 

The data used in this analysis were taken with the Belle detector [15] at the KEKB 
asymmetric energy e^e~ collider [l^. Belle is a large-solid-angle magnetic spectrometer 
that consists of a three- layer silicon vertex detector, a 50-layer central drift chamber (CDC), 
an array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement of time- 
of-fiight scintillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised 
of CsI(Tl) crystals located inside a super-conducting solenoid coil that provides a 1.5 T 
magnetic field. An iron fiux-return located outside of the coil is instrumented to detect K^^ 
mesons and to identify muons (KLM). 

The data sample consists of 140 fb~^ taken at the T(4S') resonance, or 152 x 10^ 
BB events. Another 15 fb~^ taken at 60 MeV below the resonance are used to estimate the 
jiow-BB (continuum) background. The off-resonance data is scaled by the integrated on- to 
off-resonance luminosity ratio corrected for the 1/s dependence of the qq cross-section. 
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Generic Monte Carlo samples equivalent to about three times the integrated luminosity 
are used in this analysis. Monte Carlo simulated events are generated with the evtgen pro- 
gram 171] and full detector simulation based on GEANT |18| is applied. QED bremsstrahlung 
in S — >• X^l' decays is added using the PHOTOS package [l9|. 

Hadronic events are selected based on the charged track multiplicity and the visible 
energy in the calorimeter. The selection is described in detail elsewhere |20|. We also apply 
a moderate cut on the ratio of the second to the zeroth Fox- Wolfram moment 
to reject continuum events. 



21 , i?. < 0.4, 



B. Event reconstruction 

Charged tracks are required to originate from the interaction point by applying the fol- 
lowing selections on the impact parameters in r0 and z, \dr\ < 2 cm and \dz\ < 4 cm, 
respectively. Additionally, we demand at least one associated hit in the SVD detector. For 
pion and kaon candidates, the Cherenkov light yield from ACC, the time-of-flight informa- 
tion from TOF and dE/dx from CDC are required to be consistent with the respective mass 
hypothesis. 

Neutral D meson candidates are searched for in the decays channels — >■ K~Tt^ and 

— > K'TT^TT'TT^ . We fit the charged tracks to a common vertex and reject the candi- 
date if the x^-probability is below 10^'^. The momenta of the charged tracks are re-evaluated 
at the vertex and the 4-momentum is calculated as their sum. The reconstructed mass 
is required to lie within ±3 standard deviations from m^io, where one sigma is about 4.5 MeV 
(4 MeV) for the one (three) pion mode. 

The candidate is combined with a slow pion 7r+ (appropriately charged with respect 
to the kaon candidate) to form a D*^ candidate. No impact parameter and SVD hit require- 
ments are applied for tTs. Again, a vertex fit is performed and the same vertex requirement 
is applied. The invariant mass difference between the D* and the D candidates. Am, is 
required to lie within 144 and 147 MeV. Additional continuum suppression is achieved by 
requiring a D* momentum below 2.45 GeV in the cm. frame. 

Finally, the D* candidate is combined with an oppositely charged lepton (electron or 
muon). Electron candidates are identified using the ratio of the energy detected in the ECL 
to the track momentum, the ECL shower shape, position matching between track and ECL 
cluster, the energy loss in the CDC and the response of the ACC counters. Muons are 
identified based on their penetration range and transverse scattering in the KLM detector. 
In the momentum region relevant to this analysis, charged leptons are identified with an 
efficiency of about 90% and the probability to misidentify a pion as an electron (muon) is 



0.25% (1.4%) [22|, |23|. No SVD hit requirement is made for lepton tracks. In the lab frame, 
the (transverse) momentum of the lepton is required to exceed 0.85 GeV/c (0.6 GeV/c). We 
also apply an upper lepton momentum cut at 2.4 GeV in the cm. frame to reject continuum. 
Again, a vertex fit is performed and D*~^i~ candidates are rejected if the vertex probability 
is less than 10^^. 

Figures [3] and m show the invariant mass of the candidates and the Am distributions, 
respectively. 
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FIG. 3: Invariant D'^ candidate mass distributions in the different sub-samples. All analysis cuts 
(except on the plotted variable) and | cos 9B,D*e\ < 1 are applied. 
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FIG. 4: Am distributions in the different sub-samples. All analysis cuts (except on the plotted 
variable) and | cos 6b,d*£\ < 1 are applied. 
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C. Background estimation 



Because we do not reconstruct the other B meson in the event, the B momentum is a 
priori unknown. However, in the cm. frame, one can show that the B direction hes on a 



cone around the (D*£)-axis [24 



where E"^ is half of the cm. energy and is a/ E'^ — m^. The quantities E^*^, p^,^ and 
mD*i are calculated from the reconstructed D*i system. 

This cosine is also a powerful discriminator between signal and background: Signal events 
should strictly lie in the interval (—1, 1), although - due to finite detector resolution - about 
8% of the signal is reconstructed outside this interval. The background on the other hand 
does not have this restriction. We therefore perform a fit to the cos 6 B,D*e distribution to 
determine the background normalizations from the data. 

The background contained in the selected events can be attributed to the following six 
sources: 

• continuum: any candidate reconstructed in a non-T(4S') event 

• fake lepton: the charged lepton candidate is fake; the D* candidate might be fake or 
not 

• fake D*: the D* candidate is misreconstructed; the lepton candidate is an actual 
electron or muon 

• D**: background from B z/ decays with D** D*~n or B ^ D*~Tri+u 
non-resonant 



• correlated background: background from other processes in which the D* and the 
lepton stem from the same B meson, e.g, B^ —>■ D*~t^v, r"*" ^'^vv 

• uncorrelated background: the D* and the lepton stem from different B mesons 

These background components are modeled by (i5*£)-candidates, appropriately selected 
from MC data based on generator information, except continuum which is modeled by off- 
resonance events. The shape of the fake muon background is corrected by the ratio of the pion 
fake rate in the experimental data over the same quantity in the MC simulation, as measured 
using kinematically identified pions in — >■ tt+tt" decays. The cos9b.d*i distribution in 
the data is fitted using the TFractionFitter algorithm [2^ in ROOT [26|]. The fit is done 
separately in the four sub-samples defined by the decay channel and the lepton type. 
The results are shown in Fig. [5] and Table [B 

In all fits, the continuum normalization is fixed to the on- to off-resonance luminosity 
ratio, corrected for the 1/ s dependence of the e+e~ — > qq cross-section. The normalizations 
of some other components have also been fixed to the MC expectations if their contributions 
cannot be determined reliably from the cos9b,d*i spectrum, as indicated in the table. In 
general, the normalizations obtained by the fit agree well with the MC expectations except 
for the D** component which is too abundant in the MC. 
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FIG. 5: Result of the fits to the cos 63,0*1 distributions in the different sub-samples. 



sample 


Kir, e 


Kir, fi 


KSir, e 


KSn, fj, 


signal 


(80.95 ± 1.06)% 


(80.92 ± 0.98)% 


(73.17 ± 1.71)% 


(72.22 ± 1.46)% 




(4.73 ±0.87)% 


(1.24 ±0.85)% 


(5.21 ± 1.18)% 


(2.85 ± 1.10)% 


uncorrelated 


(5.36 ±0.27)% 


(4.38 ± 0.29)% 


(5.42 ±0.58)% 


(4.17 ±0.54)% 


correlated 


(1.69 ±0.26)% 


(2.42 ± 0.28)% 


(2.04 ± 0.69)% 


(2.25 ± 0.59)% 


fake i 


0.68% (fixed) 


3.62% (fixed) 


0.72% (fixed) 


4. 04% (fixed) 


fake D* 


2.96% (fixed) 


2.91% (fixed) 


(8.78 ± 2.63)% 


(9.63 ±2.15)% 


continuum 


3.62% (fixed) 


4.51% (fixed) 


4.81% (fixed) 


4.87% (fixed) 



TABLE I: The signal and background fractions for selected events within the signal window 
\ cos Ob, D*i\ < 1- 
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FIG. 6: Reconstruction of the i? direction. Refer to the text for details. 



D. Kinematic variables 

To calculate the four kinematic variables - w, cos 6i, cos By and x ^ that characterize the 
D*~l^v decay defined in Sect. Ill At we need to determine the rest frame. The 
B direction is already known to be on a cone around the (Z}*£)-axis with opening angle 
2^s,D*^ in the cm. frame, Eq. (ITS!) . For the best guess of the B direction, we first estimate 
the cm. frame B vector by summing the momenta of the remaining particles in the event 



(^inclusive [24) choosc the direction on the cone that minimizes the difference to pfnciusive; 
as shown in Fig. [61 

To obtain pfnciusive; exclude tracks passing very far away from the interaction point 
or compatible with a multiply reconstructed track generated by a low-momentum particle 
spiraling in the central drift chamber. Unmatched clusters in the barrel region must have an 
energy greater than 50 MeV. For clusters in the forward (backward) region, the threshold 
is at 100 MeV (150 MeV). Then, we compute Pmciusive (in the lab. frame) by summing the 
3-momenta of the selected particles, 

Pinclusive = PHER + PLER " ^1 ' ^^^^ 

i 

where the indices HER and LER correspond to the colliding beams, and transform this vector 
into the cm. frame. Note that we do not make any mass assumption for the charged particles. 
The energy component of Pinciusive is defined by requiring ^in^usive to be E^^^^ = ^fsjl. 

With the B^ rest frame reconstructed in this way, the resolutions in the kinematic vari- 
ables are found to be about 0.025, 0.052, 0.047 and 6.47° for w, cosda^ cosQy and x, respec- 
tively. 



E. Fit procedure 



In the rest of this analysis, we consider only events passing the signal window requirement 
\cosQB,D*ii\ < 1- We perform a binned fit of the w, cos9i, cosOy and x distributions 
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over (almost) the entire phase space to measure the following quantities: the form factor 
normalization jF(l)|V^b| Eq- (pAj) . and the three parameters p^, -Ri(l) and -R2(l) which 
parameterize the form factor in the HQET framework Eqs. (fT2|) - f|T^ . Instead of fitting in 
four dimensions, we fit simultaneously the one- dimensional projections of w, cosOi, cosOy 
and X to have enough entries in each bin of the fit. This introduces bin-to-bin correlations 
which have to be accounted for. 

The distributions in w, cos9i, cosOy and x are divided into ten bins of equal width. The 
kinematically allowed values of w are between 1 and ~ 1.504 but we restrict the w range to 
values between 1 and 1.5. In each sub-sample, there are thus 40 bins to be used in the fit. 
In the following, we label these bins with a common index i, i = 1, . . . ,40, i.e., depending 
on the value of i, the i^^ bin might belong to the w, cosdg, cos By or x distribution. 

The predicted number of events A^*^ in the bin i is given by 

ATf = NboB{D*+ D%+)13{D°)TBori , (17) 

where Nbo is the number of mesons in the data sample and B{D*^ D^n^) is taken from 
Ref. Q. In the {K7i,e) and {K7T,fi) sub-samples B{D^) is B{D^ K-n+) [6]; in {K37i,e) 
and {K37i,n) B{D^) is Rk2,-k/k-kB{D^ — *• K~r[^), with RRSn/Kn a fifth free parameter of the 
fit. Finally, r^o is the lifetime [il, and Fj is the width obtained by integrating Eq. flTU]) 
in the kinematic variable corresponding to i from the lower to the upper bin boundary (the 
other kinematic variables are integrated over their full range). This integration is numerical 
in the case of w and analytic for the other variables. The expected number of events A^^*^^^ 
is related to A^**^ as follows 

40 

Nr = E (R^^^^K) + • (18) 

Here, is the probability that an event generated in the bin i is reconstructed and passes 
all analysis cuts, and Rij is the detector response matrix, i.e., it gives the probability that 
an event generated in the bin j is observed in the bin i. Both quantities are calculated using 
MC simulation. N^^^ is the number of expected background events, estimated as described 
in Sect. UTTCl 

Next, we calculate the variance erf of A^^''^''. We consider the following contributions: the 
uncertainty in A^**^ (poissonian); fluctuations related to A^ repetitions of the yes/no exper- 
iment with known success probability (binomial); a similar contribution to the variance 
related to Rij (multinomial); and the uncertainty in the background contribution N^^^. This 
yields the following expression for af 
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40 



i=i 



^^ - 2^ ^^3^j^3 + ^iJ^ J > + M> ^A^J > + 

^^data ^^data 



2 ^ji^ ^j) I ]yrth\2 _|_ ^ij ( ^ ^ij ) 2 / a ^th \ 2 



(19) 



+ a^(Arf'^^) . 



The first term is the poissonian uncertainty in A^*^. The second and third terms are the bi- 
nomial and multinomial uncertainties related to the finite real data size, respectively, where 
-^data (^data) the total number of decays (the number of reconstructed decays) into the 
final state under consideration {Kn or /TStt, e or ii) in the real data. The quantities ej and 
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Rij are calculated from a finite signal MC sample (A^mc and A^Mc)j corresponding uncer- 
tainties are estimated by the fourth and fifth terms. Finally, the last term is the background 
contribution a'^{N^^^), calculated as the sum of the different background component vari- 
ances. For each background component defined in Sect. IIII CI we estimate its contribution 
by linear error propagation of the scale factor and the error determined by the procedure 
described above. For continuum, we estimate the error in the on-resonance to off-resonance 
luminosity ratio to be 1.5%. 

In each sub-sample, we calculate the off-diagonal elements of the covariance matrix covjj 
as Npij — NpiPj, where pij is the relative abundance of the bin in the 2-dimensional 
histogram obtained by plotting the kinematic variables against each other, pi is the relative 
number of entries in the 1-dimensional distribution, and is the size of the sample. Co- 
variances are calculated for the signal and the different background components, and added 
with appropriate normalizations. 

The covariance matrix is inverted numerically within ROOT and, labelling the four sub- 
samples [Kir or KSn, e or /i) with the index k, the sub-sample functions are calculated, 

Xl = - NDCr/iN^'"' - , (20) 

where N°^^ is the number of events observed in bin i in the data. We sum these four functions 



and minimize the global x with MINUIT [27 



We have tested this fit procedure using generic MC data samples. All results are consistent 
with expectations and show no indication of bias. 



IV. RESULTS AND SYSTEMATIC UNCERTAINTIES 
A. Results 

After applying all selection requirements and subtracting backgrounds, 69, 345±377 signal 
events are found in the data. The preliminary result of the fit to these events is shown in 
Fig. [7] and Table HTl The statistical correlation coefficients of the five fit parameters are given 
in Table [ml 

As explained earlier, the branching fraction for the decay — > K^h^h^ti^ is floated in 
the fit to the full sample. The fit result is compatible with the recent measurement by the 



CLEO-c experiment [28|. In the sub-sample fits, Rkstt/Ktt is fixed to the value of the full 
sample fit. 



B. Systematic uncertainties 

To estimate the systematic uncertainties in the results quoted above, we consider con- 
tributions from the following sources: uncertainties in the background component nor- 
malizations, uncertainty in the MC tracking efficiency, errors in B{D*~^ —>■ D'^tc^) and 
B{D^ — > K~7r~^) j^, and uncertainties in the lifetime Q and the total number of 
B'^ mesons in the data sample. 

To estimate the uncertainty related to a given background component, we vary the nor- 
malization by ±1 standard deviation of the fit described in Sect. IIII CI The uncertainty in 
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FIG. 7: Preliminary result of the fit of the four kinematic variables in the total sample. (The 
different sub-samples are added in this plot.) The points with error bars are continuum subtracted 
on-resonance data. The histograms are the signal and the different background components. The 
color scheme is explained in Fig. [5j 

the D**lv component is inflated by a factor of two to make the amount of D**iv in each 
sub-sample consistent. The error in the continuum normalization is taken to be 1.5% as 
explained earlier. 

For the tracking uncertainty, we calculate the track finding error considering only the 
decay into K~ts:^ , as the branching fraction for K^ti^ti^ii^ is fitted from data. (A 

possible mismodeling of the tracking efficiency for this mode would be absorbed in Rkz-k/Ktt-) 
We thus have four charged tracks. Assuming 1% uncertainty for each track, except for the 
slow pion from D*^ for which we use 2%, and adding these uncertainties linearly, we obtain 
a tracking uncertainty of 5%. Given the size of this error, the uncertainty in the lepton 
identification (1—2%) can be neglected. 

The breakup of the systematic error quoted in Table [III is given in Table IIVI 



V. SUMMARY AND DISCUSSION 

We have reconstructed about 69,000 D*^i~^h'i decays using a 140 fb~^ data sample 

recorded at the T(45') resonance with the Belle detector at the KEKB accelerator. A fit to 
four kinematic variables fully characterizing this decays yields to measurements of the form 
factor normalization jF(l)|V^b| and of the parameters p^, -Ri(l) and -R2(l) that enter the 
HQET form factor parameterization of this decay. We obtain: jF(l)|\4fe| = 34.4 ± 0.2 ± 1.0, 
p2 = 1.293 ± 0.045 ± 0.029, Ri{l) = 1.495 ± 0.050 ± 0.062, i?2(l) = 0.844 ± 0.034 ± 0.019, 
and 13{B^ D*-i+ue) = (4.42 ± 0.03 ± 0.25)%. For all these numbers, the first error is the 
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odllijJic 




7T ^ jJL 


K37r^ 6 


P 


i.ozy ±u.u/z±u.ui( 


1 001 -1- n n7^ 4- n n/ifi 

i.Zzi It U.U I It u.u4D 


L.Zoo It U.Zoo It U.UOo 


I? /"I ^ 


i.400 it U.U/ / it U.U4D 


1 ^nc _i_ n nQ7 _i_ n nnn 

i.bUo it u.uo/ it u.uyy 


i.UoD it U.izo it U.U44 




n 789 -1- f) n^i^ -1- n ni 4 

u. / HI u.uotj m u.ui^ 


n QKQ _|_ n riKc; _|_ n noy 
u.ouo HI u.uoo m vj.kJz^i 


n nen -1- D 087 -1- D 097 




2.153 (fixed) 


2.153 (fixed) 


2.153 (fixed) 


^(50 ^ (%) 


4.43 lb 0.03 lb 0.25 


4.41 lb 0.03 lb 0.26 


4.42 lb 0.04 lb 0.25 


(10-3) 


34.3 lb 0.4 lb 1.0 


33.5 lb 0.4 lb 1.0 


35.6 lb 0.8 lb 1.3 




29.2/36 


37.4/36 


19.2/36 




78.2% 


40.4% 


99.0% 



sample 


KStt, ijl 


average 


full sample 


P' 






1.436 lb 0.121 ±0.062 


1.299 lb 0.045 


1.293 lb 0.045 lb 0.029 


Ri{l) 






1.643 lb 0.163 lb 0.112 


1.427 lb 0.050 


1.495 lb 0.050 lb 0.062 


^2(1) 






0.842 lb 0.105 lb 0.038 


0.844 lb 0.034 


0.844 lb 0.034 lb 0.019 








2.153 (fixed) 




2.153 lb 0.011 






(%) 


4.47 lb 0.04 lb 0.26 


4.43 lb 0.02 


4.42 lb 0.03 lb 0.25 


nmcb 


(10-3) 




35.6 lb 0.7 lb 1.3 


34.5 lb 0.2 


34.4 lb 0.2 lb 1.0 








17.9/36 




138.8/155 








99.5% 




82.0% 



TABLE II: The results of the fits to the sub-samples, their average and the fit result on the total 
sample. The first error is statistical, the second is the estimated systematic uncertainty. The 
breakup of the systematic uncertainty is given in Table lllli All numbers are preliminary. 





m\Vcb\ 


p' 


R,{1) 


«2(1) 


Rk3tt/Ktt 


mm 


1.000 


0.635 


-0.285 


-0.220 


0.011 


p' 




1.000 


0.388 


-0.870 


0.040 


i?i(i) 






1.000 


-0.511 


0.001 


^2(1) 








1.000 


0.002 


RkStt/Ktt 










1.000 



TABLE III: The statiscal correlation coefficients of the five parameters in the fit to the full sample. 



statistical and the second is the systematic uncertainty. These results are compatible with 
the recent BaBar measurements of these quantities 0]. All numbers are preliminary. 
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U.UUo 
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U.UUi 
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n 991 

U.ZZl 


u.ou 












n nsi 

U.UOi 
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U.Oi 


B{D*+ D^TT 


+) m 








0.033 
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rm m 
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/+-//00 f6] 
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0.011 


0.005 
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TABLE IV: The breakup of the systematic uncertainty in the result of the fit to the full sample. 
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